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Abstract 

Due to Girard's (sometimes called Waring's) formula the sum of the r— th power of the 
zeros of every one variable polynomial of degree N, Pn{x), can be given explicitly in terms 
of the coefficients of the monic -P/v(x) polynomial. This formula is closely related to a known 

— 1 variable generalization of Chebyshev's polynomials of the first kind, Tr^ ^\ The 
generating function of these power sums (or moments) is known to involve the logarith- 
mic derivative of the considered polynomial. This entails a simple formula for the Stieltjes 
transform of the distribution of zeros. Perron-Stieltjes inversion can be used to find this 
distribution, e.g. for N ^ co. 

Classical orthogonal polynomials are taken as examples. The results for ordinary Cheby- 
shev T]\f{x) and Un{x) polynomials are presented in detail. This will correct a statement 
about power sums of zeros of Chebyshev's T— polynomials found in the literature. For the 
various cases (Jacobi, Laguerre, Hermite) these moment generating functions provide solu- 
tions to certain Riccati equations. 



1 Introduction 

Sums of powers of the zeros of every one variable polynomial P/v(x) of degree can be 
given in terms of its coefficients as (A^ — l)-fold nested sums. If the polynomial at hand has 
definite parity, i.e. if P/v(— a;) = (— l)^P/v(x) , only [y — Ij nested sums are needed. These 



2 

„th 



facts follow from Girard's formula ||rT|: The considered r power sums are one part symmetric 
functions (polynomials) of weight r in the zeros, and due to the fundamental theorem of symmetric 
functions, e.g., ^ they can be uniquely expressed in terms of the elementary symmetric functions 
{cri}i of the zeros, i.e. (modulo signs) in terms of the coefficients of the corresponding monic 
polynomial P/v(x). Girard^s formula which solves this problem is usually, e.g., |TT| given in a 
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'semi-explicit' form where still partitions have to be performed. A simplification of this formula 
has been given in 0, [Q, [|1^. Here we shall present the explicit formula. The latter authors 
defined a certain {N — 1) variable generalization of Chebyshev^s polynomials of the first kind, 
TJ^~^)((Ti, (Tat), by the Girard- Waring expression after putting ctat, the product of the zeros, to 
one. 

The zeros of characteristic polynomials P/v(x) determine the eigenvalue spectrum of (finite-dimen- 
sional) operators. Therefore information on these zeros is welcome. In most cases these zeros are 
not known, and numerical or perturbative methods are employed. Often, the limit case — »• cx) 
is of interest and the limit density of zeros p{x) := limN^oo p(yN,x) is asked for. Here the 
power sums of zeros become important because, for given Pn{x) with only real zeros (hence for 
all orthogonal polynomials with positive moment functional), the (ordinary) generating function 
of all (positive integer) power sums, G{N,z), is related in a simple way to the Stieltjes-transfoim 
of the discrete distribution of zeros p{N, x) (the counting measure of the zeros which lives on 
the real line). This function G{N,z), which generates the generalized Chebyshev polynomials of 
the first kind, is essentially given by the logarithmic derivative of the polynomial P/v(x) under 
consideration. If it can be considered as a function of real N, and the limit N ^ oo exists, then 



one can obtain the limit distribution of zeros of Pn{,x) by Perron- Stieltjes inversion [pTSf, 



In this report which is based on the more detailed work |^ (which is referred to for proofs 
and intermediate steps) we demonstrate the above mentioned facts for the classical orthogonal 
polynomials. The generating functions G{N, z) are computed, and they are seen to satisfy certain 
i?zccatz-equations, or their equivalent linear second order differential equations. For Chebyshev's 
polynomials of both kinds simple expressions for this generating functions of the moments of 
the dicrete counting measure of the zeros result. It becomes transparent why in these cases the 
1/N expansion for the power sums or moments becomes unreliable. Finally the well-known limit 
distribution of zeros for the Jacofti-polynomials and (appropriately scaled) generalized Laguerre- 
and i/ermzte-polynomials are found from the Stieltj es-tiansform via Perron-Stieltjes-iirveision. 

As explained above, this method to find limit distributions of zeros works for all polynomials 
with real zeros whose logarithmic derivative can be considered as function of the degree number 
A^ with existing limit N ^ oo (possibly rescaled). 



2 Girard- Waring formula for the power sums of zeros of 
polynomials 

A monic polynomial of degree A^ in one real variable x is written as 

Pn{^) = E (-1)' ^k{N) x^-^ = n (x - xS^)). (1) 

fc=0 i=l 

The fundamental theorem of algebra has been used. The second equation defines the fc*'' elemen- 
tary symmetric polynomials (Jk{N) = of the A^ zeros {Xj-^''}f^ of Pn{,x) . (To{N) = 1 . 



2 



The normalized sum of the r^^ power of the zeros is 

1 ^ 

^r{N) T7 E i^^'^r (2) 

-'^ i=l 

-oo 

x"^ dpN{x) (3) 

— oo 

1 ^ 

dpN{x) := —Y,S{x-xf')dx. (4) 

i=l 

The discrete counting measure for the zeros, Pn{x), is used only for polynomials with only real 
zeros, in order to have this measure normalized to 1. Under this assumption the power sums are 
the moments of this real discrete measure. (Sometimes we shall call the power sums of zeros of 
general polynomials also moments.) 

For an arbitrary polynomial Girard^s formula, when written explicitly, looks for r G N , N = 
2, 3, ... as follows: 



m. 



W = ^ (ntY(Ei:i ^^^^)) <iv,r,{z,}f-i> nJLi^(|^^r-^ (5) 

=: i^\a^,...,ar,) , (6) 

where the N — 1 sums are ordered from the left to the right with increasing k value, and the upper 
boundaries of these sums are the greatest integer of 6i :— r/N and 

rE?-l (N+1-j) i. 
- N+l-k • * = V...iV-l. (7) 

The combinatorial factor which appears in eq. (5) is 

• (r-Ef.l'(iV + l-j)i,)! ■ 

For r = one has mo(A^) = 1. t^^^ai, ciAf), which is a polynomial, has been introduced because 
it is related to an iV — 1 variable generalization of Chebyshev's polynomial of the first kind T^^x). 
This connection is given after the rescaling 

r(^-^)(.i,...,s^_i) =iW(<7i,...,a^)/Kr/^ , (9) 

with 

Sk ■^(Jkl{<JNf''' , A;=l,...,7V-l . (10) 
The ordinary Chebyshev polynomials are now Tj.{x) — T^^\2x). 
Corolleiry: Moments for Pn(— x) = (—1)^ Pn(x) polynomials 

i) mo{N) = 1 , AT e No a) m2i+i{N) =0 , / e Nq , iV e N , (11) 

in) m2i{N) = 
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= ^LfJ(-iy^!^^^^(^2,a4,...,a,L^j) (13) 
= FLfJ(-iy ^Lfj)'^^^^ T/L^J-^\ei,...,eL.j_,) (14) 

for /e No, AT = 4, 5,..., 

[—1 — 1 (l~J) 

where Bk, resp. < [yj,^ > , resp. ^ , is obtained from bk of (7), resp. 

< N,r, > of (8), resp. t^^^ of (6), by replacing — > [yj and r ^ I. The variables of 

Ti ^ , which is also obtained by the same replacement from T^^"^) given in (9), are 

AT 

In addition: 17121(1) = , ma^S) = (-^^2)' , m2/(3) = §(-^2)^ 
Example: 

Consider Chehyshev^s polynomials of the first kind, Tn{x), which are classical orthogonal poly- 
nomials of the Jacobi type f^. The zeros are known, and the non- vanishing moments are 

m^\N) = Z^iii (cos (^I^tt))^' . The Corollary gives an alternative expression for this sum. 

For/eN,MGN : m^P(2M) = T/*^"^^(ei(2M), eM-i(2M))/2(2-i?)' , (16) 

with efc(2M)^^2-^(2*^,:^,-i) , fc = l,...,M-l, (17) 
m^P(2M + l) =^ ((H^^)'Tf^-^)(ei(2M+l),...,eM-i(2M+l)), (18) 

with e,i2M+l)^^-^^-^{^l~'), A: = 1,...,M-1. (19) 

andmf^(A^) =1 , m^^\l) = 0. 
Some instances for / G N are: 

m(P(2) =(!)' , m(P(4) = TKv^)/(2v^)' , 

m^P(6) =T/'^(3-2l,9-2-i)/2i' , (20) 
mS\3) =3'-V22,-i , m(P(5) = (v^)'-^ TKv^/2)/22{'-i) , 

^^^(7) =f (4yT/^^(7i,2-7i) . (21) 



From this one finds for up to 7: 



N=l 
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These values are in conflict with a statement found in There it is claimed that the moments 
(22) are given by the highest values, independently of A^, which is false. Only the values for all 
N > ^ + 1 are correctly given in [Q. 

The example of Chehyshev's polynomials of the second type is also treated in 0. 



3 Moment generating functions 



The generating function of Pr{N) , the sum of the r*^ power of variables {xi}^ , is given by the 
logarithmic derivative of the generating function of the elementary symmetric functions {ck} of 
these variables, see e.g. p. With 



N 



N 



E{z) ^ E^{z) := n (1 - ^) = E i-^r^r z 

i=l r=0 

which is also given by E{z) = z^ Pj^{l / z) , one finds from 

{lnE{z))' =E\z)/E{z 
and Gimrd' s iomuXa, (5), (6) 

oo 



(To = 1 



-| oo 
Z r=l 



(23) 



(24) 



G{N,z) ■.= Y.mr{N)z' 

r=0 



1 ^l^ivW 



N Pn{x) 



x=l/z 



(25) 



As an application we computed the moment generating functions G{N, z) for the classical orthog- 
onal polynomials: Jacobi P, 



(",/3), 
N ' 



X) 



generalized Laguerre L^^\x) and Hermite Hn{x) . The 
derivative rules for these polynomials (c/. |^) have been used. We quote the results for Chehy- 
shev's polynomials of both kinds separately because they can be written in terms of an elementary 
function. For details see The hypergeometric function 2-^1(0, b, c; x) and its confiuent version 
F{a,c;x) = iFi{a,c;x) appear in table 1. 

The (A^ + 1)— term recursion relation for the 4^^— polynomials which can be derived from eq.(25), 
is: 



N 



for n 



N,N + 1,...,: 

j=0 







(26) 



where the input quantities {tj }q ~ are either computed from Girard's formula or from the 
generating function G{N, z). 

We also comment on the generating functions G'< ( A^, z) for the normalized sums of negative 
powers of the zeros: m_r{N) := j^Y^^i [xf^^)'"^ , r G N. Here we restrict the discussion to 

polynomials which have no vanishing zeros. Withii^^(2;) := YlfLi{l—z/x\^'^) = {—1)^Pn{z)/(Tn 
one finds 

z P'^iz) 



G4N,z) := J2 m-r{N) z' 

r=l 



-j^iln E<{z))' 



NPn{z) 



G{N, 11 z) . (27) 
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Beware of expanding the identity G^{N,z) + G{N,l/z) = in 2 and 1/z (and concluding 
falsely that all power sums have to vanish). The sums of negative powers of the Pn{x) zeros can 
then be found by differentiation. The connection between Fad di Bruno^s formula and the Bell- 
polynomials |]13| and the identification of the generalized Chebyshev polynomials t^^^ as special 
5e//-polynomials leads, with the abbreviation g{x) = Pn{x) , for m G N first to 

^ In g{x) = 

-{m-l)\ N t^^\-g^^\x)/{l\g{x)), +g^^^{x)/{2\g{x)), ...,{-l)"'g^"'\x)/{m\g{x))) , (28) 

with (yf^^) = ^g{x) . After putting x = 0, using {—lY g^^\Q) / {k\g{Q)) = dN-k/f^N , one finally 
finds 

{ti^\crN-i/crN, crN-2/crN, crN-r/crN,0, ...,0) ioi < r < N 
(29) 
4^)((TAr_i/aAr, (TN~2/crN, (<To = l)/a7v) for r > N , 

where in the upper alternative N — r zeros have to be inserted. 

Table 1 : Generating functions for power sums of zeros of classical orthogonal 

polynomials 



Pn(x) 


G(N,z) 


p(",/3) 


(^{a,0)(AT ^\ N+a+l3+ll2Fi(-N+l, N+a+13+2, a+2; {z-l)/2z) 
l^^'^J 2{a+l) z 2Fx{-N, N+a+P+1, a+l; (z-l)/2z) 


Tn 


G(^)(iV,z) - -^tanh{N In^^) 


Un 


G(^)(iV,z) - -^[{l + j^)tanh({N + l)ln^^)-j^^] 


r(") 




Hn 


riH)(AT .^ _ in-V, h iA^)/r(2^) -1 i^(2^. I; V^')/rC-^) 



One can verify (c/. 0) that for orthogonal polynomials the generating functions G{N,z) are 
special solutions of certain iiiccati-equations which we list in the following table 2. 
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Table 2 : Riccati-equations satisfied by moment generating functions of classical 

orthogonal polynomials 



■p 


Riccati-eq. for G(N, z) for every N 




+(1 


z) = + (a - + a + /5 + 1) ;^ G^^'W (iV, z) 


Tn 




) ^^iG^^\N,z) ^{\-z^){G^^)[^N,z)f -1 


Un 






r(«) 
J^N 




= iV^(G'(")(Ar,2))2 +(q;2_1) ^("^(A^,^) +1 


r J") 


z'^-§:Q(-HN. z) 


= ^(Q(«)(iV,2))^ + (a^/iV- 1) g(°)(iV,^) +1 


Hn 


f£Gi^\N,z) 


^^7V(G'(^)(7V,^))2 _(i + |)G'W(7v,^) +1 


Hsn 




= 4(g(^)(iV,^)f -(i + ^)gW(iv,z) +1 



We have also listed the i2iccatz-equations satisfied by the moment generating functions of scaled 
generalized Laguerre- and i^^erm^^e-polynomials. The definitions are: 

Ls^^\x) L^S\Nx) , Q^^\N,z) := G^^\N, z/N) , (30) 

Hsn{x) := Hn{\/Nx) , Q^^\N,z) := G^"\N , z / y/N) . (31) 

The associated second order linear differential eqs. are listed in table 3. G{N, z) , 
resp. g(")(iV,;2) := G^") (iV, ^/iV) and Q'^"\N,z) := G^"\N,z/VN), are replaced by 
~'k ^(/n U{N, z)), with correspondingly labelled U{N, z). 
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Table 3: Second order linear differential equation satisfied by W(N, z) for given N for 

(scaled) classical orthogonal polynomials 



Pn(x) 


2°<^ order LDE for W(N, z) for every N 




(^2(1 -z^)dL- z{2z^ + (a - + a + - N{N + 1 + a + f3)) U'^^^^^N, = 


Tn 


(^2(1-^2)^ 2^2)1 -iv2)w(^)(iV,z) =0 


Un 


/ O/-. 0\/-l2 ^0\^^ -7t-r/-7t-r ^\\-»i/rr\/-7t-r \ „ 

(2;2(l-2;2)^ _2;(l + 22;2)£ - Ar(iV + 2)jW(^)(Ar,2;) =0 




(z3 ^ +;2(Ar+(l-«)z) ^ + Ar2) w(")(iV,z) =0 


Hsn 


(;24 ^ +2z(z2 + Ar)^ +2Ar2) wW(Ar,2) =0 



Tfie above given Riccati-eqyiaiioii's can also be found starting with Case's eqs. which have been 
derived for the sums of powers of (simple) zeros of polynomials satisfying certain differential 
equations [||. Table 4 shows these eqs. for the classical orthogonal polynomials. In the generalized 
Laguerre and the Hermite case scaled moments are used in accordance with the above defined 
scaled polynomials Ls^^\x) and Hsn{x) : 

gf)(iV) =m(^)(iV)/iV^/2 ^ q(»\N) =mi''\N)/N' . (32) 
The input moments are: 

Jacobz: moiN) = 1 , mi{N) := {/3 - a)/{2N + a + f3) . 
FotTn : a = (3 = -1/2 and for Un a = (3 = +1/2 . 
generalized Laguerre: qQ°'\N) = 1 . 
Hermite: q^"\N) = 1 , q["\N) = . 

From the Riccati-eqs. which the generating functions of the power sums satisfy, one infers the 
N ^ oo limits (the leading term of the 1/A^ expansion) which are given, together with the 
(scaled) moments, in table 5. Gq{z) := limAr^oo G^°''^\N, z) in the Jacobi-case. In the general- 
ized Laguerre-case Qo{z) := lim^r^oo Q^'^^N, z) , and in the Hermite-ca.se Q\^\z) := Q^^^{N, z) . 
c{z) is the generating function of Catalan^s numbers {C^}. Universality is manifest due to the in- 
dependence of the parameters a,P, resp. a, in the Jacobi-, resp. generalized Laguerre-case. These 
asymptotic moments have also been considered in |jl2[ ( although the appearance of Catalan^s 
numbers has not been mentioned). 
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Table 4: Case's equations for moments of zeros of classical orthogonal polynomials 



Pn(x) 


Ceise-eqs. 


for {mr(N)}, resp. {qj.(N)}, for given N 


p(",/3) 


{a + f3 + 2N-r) mi-^f (AT) =(/?-«) m(°'^)(Ar) - r m^^'f (TV) + 
+ NElZlimi'-jtsiN) rr^'^KN) -m^f (iV) m£;f (iV)) 


Tn 


(2A^- 


- (r + 1)) mJ+\(A^) = -r mS(A^) + 


Un 


{2N- 


- (r - 1)) mJ'^U^) = ^S(^) + 
:m£L.(A^) mP{N) - mf_l{N) m%\{N)) 


r J") 


dim 


= ^ d^KN) +E:=o ll-siN) g(")(iV) 


Hsn 







Table 5: (Scaled) moment generating functions in the limit N — > oo 



Pn(x) 


Go(z) , resp. Qo(z) 


limit moments 


p(«,/3) 


G'o(^) =l/Vl-z^ 






Qo{z) =^(l-Vl-4z) =:c(z) 


go,. =a (t'-)/(r + l) 


Hsn 


g(^H^) =c(^V2) =^(l-Vl-2^2) 


%,2l —^l/^ 
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It is clear from the explicit expressions for the generating functions G{N, z) in both Cheby- 
shev instances that a expansion is doomed to fail. Such an expansion produces for the 
Tjv— polynomials the asymptotic moment generating function Go{z) given for the Jacobi-ca.se in 
table 5. This happens to be a solution of the Tn — Riccati-eq. shown in table 2 for all A^. However, 
it is not the correct solution G^'^\N, z) given in table 1. This explains the mistake which occured in 
Q. It is possible to recover G^^\N, z) from the sj^ecm/ solution Gq{z) of the T/v— i?«ccat2-equation. 
In the Chebyshev Un case the sum of the first two 1/A^-expansion coefficients furnishes a special 
solution of the [/at— Riccati-eq. (see 0), which is also not the correct solution G^'^\N^z) given 
in table 1. 



4 Distribution of zeros and Perron-Stieltjes inversion 

For polynomials with only real zeros the (ordinary) generating function G{N, z) 
of the sums of powers of zeros (or moments) is related to the Stieltjes-iiansioim. X{N, z) of the 
discrete real measure (4) in a simple way. The Stieltj es-transfoTm of pn{x) is defined by 

/+00 \ 
dpN^x) , z ^ supp{dpN) ■ (33) 
-oo Z — X 

Because (in the sense of formal power series) one has with the counting measure (4) 

oo 1^1 r°° 1 
G{N, z) mr{N) =YrT. HvY = / 1 dpr,{x) , (34) 

r=0 i=ll- Z X- ' -^-oo i- - Z X 

one finds the simple relation 

X(N,z) = - G(N,-) . (35) 

z z 

It is possible to recover the discrete measure from X{N, z), i.e G{N, z), with the aid of the Perron- 
Stieltjes inversion formula |T5| , P , P] , P 

pNit2) - pNih) = lim / dt ImX(N,t + iri) , (36) 

vr v^+o Jti 

with pAr(tfc) := \{pN{tk + 0) + pNitk — 0)) for /c = 1, 2. This inversion is valid for real measures 
of bounded variation, provided Xj^ = X]^{z) ioi z & and X^v is analytic off the real axis. As 
long as is finite, Xj^^^z) has poles only at the zeros of P/v. The inversion formula then yields 
for pjv jumps at these real zeros ajj^"* with height m^^^ /N where mj counts the multiplicity of the 
zero x[^\ ((36) shows that poles in X lead to jumps in p. The converse statement is obvious from 
(33).) In this way one recovers the purely discontinuous measure (4). 

The inversion formula becomes especially interesting if the limit X{z) := limAr_^oo ^(A^, ^) exists. 
Then one can find the limit measure or distribution of zeros. In our examples of classical (scaled) 
orthogonal polynomials the limit moment generating functions exist (see table 5), and hence one 
can find by Perron-Stieltjes inversion ofX{z) := i Go(-) (resp. i Qo{-)) the distribution of zeros 
in the limit —>■ oo. These are listed in table 6. In all cases the conventional definition of the 
square-root in C has been used. These limit distributions (with other normalization) have been 
considered in Usl. 
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Table 6: Limit distributions of zeros of (scaled) orthogonal polynomials 



Pn(x) 


p{x) = limAr^ooPAf(a;) 




pia,l^)(^x) — - 3 , for \x\ < 1 and otherwise 


r 


p(Ls("))j^^-j _ J-^^ — 1 , for < X < 4 , and otherwise 


Hsn 


p(^*)(x) = — , for \x\ < \pl and otherwise 



5 Conclusion 

Formula (25) which gives the (ordinary) generating function GiN ^ z) for the sums of (posi- 
tive) powers of zeros of any one-variable polynomial (1) can be simplified whenever the (logarith- 
mic) derivative rule for the polynomial under consideration is known. If the polynomial satisfies a 
(second) order differential equation (like in the assumption of the work 0) one expects G{N, z) to 
satisfy a Riccati-tjpe differential equation. For polynomials with only real zeros {e.g. orthogonal 
polynomials with positive moment functionals) one can use Perron- Stieltjes inversion to find the 
distribution of zeros in the limit N ^ oo, provided limAr^oo G{N, l/z)/z exists. 
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